Nonlinear Effects in the TGB^ Phase 
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We study the nonlinear interactions in the TGBa phase by using a rotationally invariant elastic 
free energy. By deforming a single grain boundary so that the smectic layers undergo their rotation 
within a finite interval, we construct a consistent three-dimensional structure. With this structure 
we study the energetics and predict the ratio between the intragrain and intergrain defect spacing, 
and compare our results with those from linear elasticity and experiment. 
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The resolution of frustration is a central theme in con- 
densed matter physics. Competing terms in the free en- 
ergy can favor incompatible spatial organizations, the 
resolution of which often leads to a rich phase behav- 
ior and intricate spatial patterns. The TGB phase of 
chiral smectics is an ideal example of the resolution of 
frustration |lj . While the smectic part of the free energy 
favors a lamellar structure, the chiral nematic part favors 
a uniform twist of the nematic director. These two struc- 
tures are incompatible and so the layered structure must 
be riddled with defects to accommodate the twist. The 
TGB phase balances the competing interactions by form- 
ing a lattice of screw dislocations arranged in twist grain 
boundaries. In earlier work we investigated the geom- 
etry of the TGB phase within a harmonic free energy in 
which topological defects interact via screened exponen- 
tial potentials. However, it was argued in || and B that 
defects in the same grain boundary interact via power- 
law potentials when those nonlinearities required by rota- 
tional invariance were included in the free energy. In this 
letter we investigate the rotationally- invariant, nonlin- 
ear energetics of the full three-dimensional TGB^ phase, 
where the smectic blocks between the grain boundaries 
are smectic- A. We find that intergrain and intragrain 
interactions of dislocations are both power-law, and we 
compute the aspect ratio between the intragrain spacing 
Id and the intergrain spacing 4. Our computed value is 
of the same order of magnitude as that computed from 
the linear theory and measured by experiment. We will 
comment on the discrepancy and its source in the con- 
clusion. 

In an earlier paper ||, we considered the screw dislo- 
cation lattice in the TGB^ phase near the transition to 
the nematic phase. In this regime the dislocation den- 
sity is low, so the interaction between dislocation cores 
is negligible. This implies that the lattice structure is 
completely determined by the elastic energy cost of the 
smectic distortions created by the screw dislocations. We 
demonstrated in || that these distortion fields can be 
consistently treated within the harmonic approximation 
to the elastic free energy. In this approximation, we ex- 
press the elastic free energy in terms of the (Eulerian) 
layer displacement field u and the deviation <5n of the 



nematic director from the z-axis: 
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+Ki(V± • <5n) 2 + K 2 (V± x Sn) 2 + K 3 (d z 5n) 2 

where B and D are elastic moduli and K\ , K% and K3 are 
the splay, twist and bend Frank constants, respectively. 
This free energy is analogous to that of a type II super- 
conductor and the topological defects in a smectic are 
likewise analogous to Abrikosov vortices. Employing (|l|) 
the interaction energy between two parallel screw defects 
is 
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where /Co is the modified Bessel function of order zero, 
d is the equilibrium layer spacing of the smectic phase, 
and A = y/K^jD. 

A twist grain boundary is composed of an array of 
equidistant parallel dislocations, separated by a distance 
Id, which force the smectic layers to rotate. On average, 
the smectic layers drag the nematic director along and 
allow it to twist. As with screw defects and vortices, 
the TGB phase is the analog of the Abrikosov vortex 
lattice phase in superconductors [Q, |) , although the lat- 
tice structure of defects in the TGB phase is significantly 
more complicated due to the relative rotation of disloca- 
tions in different grain boundaries. We showed in (|] that 
the linearized interaction energy of two grain boundaries 
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Thus the linear theory predicts an exponential interac- 
tion between screw dislocations both within the grain 
boundary and between different grain boundaries. It is 
simple to check that the elastic interactions of disloca- 
tions within a grain boundary are balanced by the in- 
terboundary interactions, which leads to a stable ratio 
R = ld/h ~ 0.95 Q, within the experimental range of 
0.74-1.08 § §. 
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Since both smectic and cholesteric ordering are spon- 
taneously broken symmetries, it is essential that our free 
energy be invariant with respect to rotations. The har- 
monic free energy is only invariant under infinitesimal ro- 
tations; the extra terms needed to render the free energy 
(|l|) invariant under finite rotations lead to nonlinearities 
that dramatically change the energetics of even an iso- 
lated screw dislocation j|. In order to investigate these 
effects, we work with a rotationally invariant theory that 
depends on the only Goldstone mode in the system, the 
layer displacement field u jjj: 

F = \j d?x {B(v z + v 2 /2) 2 + Jf (V- v) 2 } , (4) 

where v = Vu is the layer tilt field, which, on average, 
follows the deviation of the nematic director. The free 
energy (m) includes director modes and the interaction 
energy accounts for those modes outside the core. In 
the nonlinear theory that we are now considering, these 
director modes are absent. We believe that the essential 
physics is unchanged but that the energetic details will 
require incorporation of these modes |f 2| . Note that we 
have not used the mean curvature H = V [(z+v)/|£+v|], 
but rather V-v in the bending term. Both expressions are 
rotationally invariant and though the former has greater 
geometric significance || , the latter is significantly easier 
to manage. In the limit of low-angle grain boundaries, 
the two expressions concur. Moreover, for the energy 
minimizing structures we will consider here, the mean 
curvature remains finite everywhere as long as V • v is 
finite. 

We construct a grain boundary as a linear superposi- 
tion of parallel equidistant screw dislocations (LSD) [Q, ||] 
and then use it to construct a TGB^ phase by com- 
pressing each of these individual boundaries along the 
pitch axis to force their rotation into a finite interval. 
We may then piece together a periodic lattice of these 
grain boundaries through translation and rotation in or- 
der to construct a full, three-dimensional TGB structure. 
We will calculate the free energy of this structure using 
(^J) to find a minimum energy solution and consequently 
the geometric parameter R = ld/h for the system. As 
noted in Q| , an important difference between the nonlin- 
ear and the harmonic theories is that in the former the 
interaction between the dislocation cores may not be ne- 
glected: the size and shape of the core region is another 
variational parameter. 

We start with the layer distortion field for a linear su- 
perposition of screw dislocations (LSD) in a smectic liq- 
uid crystal. Through deformation this will become the 
displacement field in a single smectic- A block. Choosing 
the average layer normal to point along the z-axis, and 
expressing the displacement u and coordinates x and y 
in units of ld/2ir, we have M 
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tan(y/2) 
l)z. 



The layer rotation angle a across the grain boundary is 
set by the topology of the grain boundary and, there- 
fore, does not depend on specific details of the grain 
boundary model. It is determined by the layer spacing d 
and the dislocation spacing Id through sin(2-) = d/2ld- 
The first term in (|5|) describes an array of parallel screw 
dislocations aligned along the z-direction at i = 0, 
y = 0, ±ld, ±2ld, ■ ■ ■ , while the second term ensures that 
the nonlinear compression energy completely vanishes far 
away from the dislocation array. 

To construct the TGByt phase out of individual twist 
grain boundaries, we confine the layer rotation of the 
regular LSD to a finite interval, (— /, I) so that the grain 
boundaries are spaced 21 apart in units of Id /lit, so 
21 = Y~lb — 7f ■ To do this, we replace the coordinate x in 
(||) with an odd deformation function, fi(x) which mono- 
tonically maps the finite interval (—1,1) to (—00,00). So 
that we recover the single grain limit, fi(x) = x for large 
I. The energetics of the confined LSD will impose ad- 
ditional constraints on fi(x). These constraints and the 
optimal choice of fi(x) will be the focus of our calcula- 
tion. Since the energy depends on v we consider: 
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V 2 / cosh fi (x) — cos y 
v z = cos (!)-!. 
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The compression elastic free energy per dislocation per 
unit length in the z-direction can be found using (§) and 
the expressions ([|||): 

^ = 2W (|) (^) 3 Wtf.a;*], (9) 

where J compr is a dimensionless functional that depends 
on the deformation function, a short distance cutoff a(x) 
required to excise the core from the integration domain 
|Q, 0], and the lattice aspect ratio R. It is interesting 
that this energy is essentially nonlinear since it depends 
on the fourth power of the displacement field u H . 

On the other hand, the curvature contribution to the 
free energy per dislocation need not have a short-distance 
cutoff: as we will show, by judicious choice of the defor- 
mation fi(x), the curvature is everywhere finite. Inte- 
grating along the grain boundary direction (y) we find: 
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Since a single helicoidal screw dislocation has no curva- 
ture, we need not excise the core region in the curvature 
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energy as long as we maintain the screw-like character of 
the defect. Near x = 0, fi{x) » is small and so we may 
expand the integrand X C urv : 
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It is evident from the above expression that the curvature 
energy diverges near x = unless //(0) = 1 and //'(0) = 
0. Thus we have two additional constraints on fi(x). To 
finish our analysis we must consider the behavior of the 
energy halfway between grain boundaries at x = ±Z. It 
is there also that we must take care to match the layer 
normals in successive blocks so that the nematic director 
is well defined everywhere. 

Near x = ±Z, the compression energy is finite since v 
is well behaved. However, the curvature energy is more 
problematic. Since fi(x) — > oo as x — > I, we can expand 
the integrand in powers of e~' x : 
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Since }i{x) diverges at s = I, we might consider set- 
ting the largest term to zero, i.e. the n = 1 term in 
(|l2|). Implementing the boundary conditions /;(0) = 
and fi(l) = oo, we find: 



fiix) = sgn(a;)log 



sinh I 



sinh(Z — \x\) 



(13) 



It is straightforward to check that as x — > I, the only 
nonvanishing term in (||) is the n = 2 term and that at 
x = I 
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We can embellish (|l3|) with a power series around |a;| = I 
so that the first and second derivatives satisfy the re- 
quired conditions at the origin: 



fi(x)= sgn(ar)log 



sinh I 



sinh(Z — \x\) 
+sgn(a;)[a + a 2 {l - \x\f + a 4 (l ~ \x\)% 

where ao , ai and 04 may be straightforwardly calculated 
JTo| ]. These extra terms do not spoil the divergent be- 
havior at x = I and so it would appear that we have 
found a low curvature deformation of the original grain 
boundary. However, our goal is to assemble individual 
blocks into the TGB structure. This is only possible if 
the smectic layers match at the midplane. Computing 
the layer tilt at the midplane, we find that v x has a small 
periodic component, while v y attains a constant value 
corresponding to the rotated layers: 

% / a \ sinw , „ 

v x (l,y) = -2am (-)^, (16) 

v y (l,y) = sin (|) . (17) 



Since v x reflects the periodicity of a single dislocation 
array and in adjacent blocks the dislocation arrays are 
rotated with respect to each other, the nonvanishing of 
v x at the midplane makes a perfect fit impossible. This 
might be a consequence of the limited class of transfor- 
mations considered in our model that allow the smectic 
layers to relax only in the direction of the pitch axis. A 
more general model which allows the defects in the layers 
to ripple in harmony with neighboring grain boundaries 
might allow for a nontrivial value of v x at x — ±Z. In- 
deed, rippling defects have been considered to explain 
commensurate twist grain boundary structures u% . 
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FIG. 1: Dependence of I CUIV (R;m) on the deformation pa- 
rameter m. 

The smectic layers could be matched perfectly if fi 
were modified to diverge faster than the form in (113). 
Consider a one-parameter family of deformations based 
on ©: 



( sinh(l/m 
fl(x;m) = m sgn(a;)log | ^T^77j 
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As before, these deformations can be altered so that they 
satisfy all the necessary boundary conditions. As \x\ — > I, 



(15) both f['(x; m) and [f[(x; to)] diverge as (I- 



) , while 

exp{— 2fi(x; to)} ~ (I — |x|) 2m . Examining I curv , we see 
that as long as to > 2, the curvature remains finite in 
the whole region. In addition, v x <~ (/ — |x|) m_1 , and 
so for v x to vanish at x ± I, we must have m > 1. By 
numerically evaluating icurv, we find that the curvature 
energy is an increasing function of to, as shown in Figure 
0. Thus we choose to = 2 in our calculation of the aspect 
ratio R. We note that the isolated m = 1 deformation 
has a still lower energy, though it is, unfortunately, not 
allowed because of geometry. In future work we will re- 
consider the to = 1 deformation by allowing for more 
general variations of the dislocation lattice. The equilib- 
rium lattice configuration is determined by minimizing 
the total free energy density with respect to the dislo- 
cation spacing, core size, and parameters characterizing 
the core shape. The total energy density includes the 
elastic energy we have discussed, an energetic cost from 
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ergetics [H. In computing I c , 



the core size has to 



be reexpressed in units of Zrf/27r, so that it does not vary 
as Id is varied. In units of ld/2-K, the core size is a/y/R. 
The entire minimization procedure can be formulated in 
terms of the dislocation density 1/c 2 , the aspect ratio of 
the dislocation lattice R, and the core size a. Recalling 
that sin (2-) = d/2ld and assuming K/B = d 2 , we have 
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FIG. 2: Dependence of the optimal lattice aspect ratio R and 
core size a on the control parameter n 2 /e for m—2. 



the dislocation cores and the crucial chiral contribution 
which favors twist. We assume that the core energy is 
proportional to the cross-sectional area A of the core re- 
gion, with an energy per unit area E. The chiral energy 
gain per dislocation H is independent of the details of 
the dislocation arrangement, and so for given values of E 
and H we must minimize 
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with respect to Id, h, and A. Note that the integrals ap- 
pearing in the compression and layer curvature terms do 
not depend explicitly on the layer rotation angle a. It is 
convenient to separate the effects of varying the disloca- 
tion density \/{ldh) = 1/c 2 and the lattice aspect ratio 
R = ld/h- In terms of R and c, Id = c\J~R and lb = cj \f~R~. 
For simplicity, we assume that dislocation core regions 
are square with sides a in units of c/2tt. Our numeri- 
cal investigation of the I C ompr[ fi(x), a(x);R] suggest that 
the core, though essential, does not greatly alter the en- 



where e = (32tt 2 / Bd A )E and T) = (32tt 2 / Bd 4 )H . When 
we minimize /tot with respect to c, a and R, we find that 
the equations for the optimum values of R and a only 
depend on the combination rj 2 /e. The optimal values of 
R and a computed for m = 2 are given in Figure |^. 

We find that the variation in the lattice aspect ratio 
mainly occurs for large core energies and then asymp- 
totes rapidly to R ~ 0.46 (see Figure^). In comparison 
with the linearized theory which predicts that Id ~ h, the 
model studied here predicts that the repulsion between 
grain boundaries is stronger than the repulsion between 
the defects in the same grain boundary and hence h is 
roughly twice as large as Id- This may be a consequence 
of our locking the director to the layer normal near the 
defect cores Q. Indeed, from Figure |we see that the op- 
timal cores are rather large, about 1/6 the defect spacing. 
In future work we will reintroduce the director modes in 
a rotationally invariant fashion. This will certainly lower 
the overall energy and should allow the cores to shrink. 
As a final aside, we note that were we to choose m = 1 in 
( p~8| ) and simply ignore the director mismatch, we would 
find R « 0.85. That solution will be studied in further 
work @. 
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